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■ Abstract. The present paper has a twofold contribution: first, we intro- 

duce a new concept of Hardy spaces on a multidimensional complexified 
annular domain which is closely related to the annulus of the Klein-Dirac 
quadric important in Conformal Quantum Field Theory. Secondly, for 
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harmonic Gaufi-Jacobi cubature formulas, which have been introduced 
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1. Introduction 

In the classical one-dimensional theory of quadrature formulas there are dif- 
ferent approaches to the estimation of the quadrature formulas of Gauss type. 
The first approach established already by A. Markov estimates the error of 
a quadrature formula for differentiable functions in C N (I) defined on the 
interval I by means of its N— th derivative, cf. Krylov [31], chapter 7.1, or 
Davis |10t p. 344] . However this approach is usually not very practical beyond 
derivatives of order five, see [52]. The second approach estimates the error 
for certain classes of functions / which are analytic on some open set D in C 
containing the interval /, cf. [TT], chapter 4.6, see also [3T], chapter 12.2. 



The first-named author thanks the Alexander von Humboldt Foundation and Project DO- 
02-275 with Bulgarian NSF. 
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Approximation of integrals in the multivariate case is a much more 
difficult task. In Numerical Analysis, instead of quadrature formula the notion 
of cubature formula is often used, see [37], [H]> [3H] and the recent survey 
[9]. In contrast to the univariate case there is no satisfactory error analysis 
available in the multivariate case, cf. [51], [6], and part 4 in the last Russian 
edition of the classical monograph [31 . In passing we mention that the area 
of quadrature formulas for harmonic functions in which Werner Haufimann 
has worked, are an interesting multidimensional alternative and we refer to 
[2], [H] on this topic. 

1.1. Preliminaries 

In |23j and |21] we have introduced a new multivariate cubature formulae 
Cn if) depending on a parameter N g N which approximates the integral 

J f(x)dp(x) (1.1) 

A a , b 

for continuous functions / : A a _b — > C defined on the annular region 

A a ,b = {x e R d : a < \x\ < b} , (1.2) 

where |x| = y/x\ + .... + x 2 d is the euclidean norm of x — (xi, ...,Xd) & R d . 
In [33], [3T], we have considered a special class of signed measures /i with 
support in A a ^, the so-called pseudo-positive measures (see definition below) 
however this formula is easy to extend to almost all signed measure. The 
exact definition of the cubature formula CV (/) will be explained in Section 
[71 One aim of this paper is to provide an error analysis for a class of functions 
on the annular region A a ^ which exhibit a certain type of analytical behavior. 

In order to give the reader the needed mathematical background we 
have to recall some terminology. Let 

S^ 1 := {x e M d : |x| = 1} 

be the unit sphere endowed with the rotation invariant measure d9. We shall 
write x € M. d in spherical coordinates x = r9 with 9 g Let T~Lk (R d ) be 

the set of all harmonic homogeneous complex- valued polynomials of degree 
k. Then / G H k (K d ) is called a solid harmonic and the restriction of / to 
S d_1 a spherical harmonic of degree k and we set 

a k :=dim?4 (M d ) , (1.3) 

see [50], [45], pQ, [20] for details. Throughout the paper we shall assume 

r M :R d ^R,^ = l,...,a fc , (1.4) 

is an orthonormal basis of Hk (R d ) with respect to the scalar product 

(/,0>s*-i == / f(0)W)dB. 

We shall often use the trivial identity Y k ,i (x) = r k Yki (0) for x — rO. The 
class of pseudo-positive measures used for our cubature formula Cn (/) is 
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now denned in the following way: a signed measure /i with support in A a ^ C 
R d is pseudo-positive with respect to the orthonormal basis Y k ,i,£ — 1, ...,<2fc, 
fc G No if the inequality 

/ /i(|x|)y M (a;)d/i(x)>0 (1.5) 

JR d 

holds for every non-negative continuous function h : [a, b] — > [0,oo) and for 
all k € No, £ = 1,2, ...,afe. Let us note that every signed measure d/i with 
bounded variation may be represented (non-uniquely) as a difference of two 
pseudo-positive measures. We refer to [33] for instructive examples of pseudo- 
positive measures. 

The cubature formula Cn (/) approximating the integral (jl.ip is based 
on the Laplace-Fourier series of the continuous function / : A a ^ — > C, 
defined by the formal expansion 

oo a k 
fc=0 i=l 

where the Laplace-Fourier coefficient ff./ (r) is defined by 

/ W (r) = / f(r9)Y k ,t(6)d6 (1.7) 

for any positive real number r with a < r < b and is defined in (jl.3l) . 
There is a strong interplay between algebraic and analytic properties of the 
function / and those of the Laplace-Fourier coefficients f k ,i- For example, if 
/ (x) is a polynomial in the variable x = (xi, xj) then the Laplace-Fourier 
coefficient fk/ is of the form fk y e {r) = r k pk,e (r 2 ) where p k ,e is a univariate 
polynomial, see e.g. in |50j or |47j . Hence, the Laplace- Fourier series (jl.6p of 
a polynomial / (x) is equal to 

deg/ a k 

f(x)=J2J^p k ,i(M 2 )Y k , t (x) (1.8) 

fe=0 £=1 

where deg / is the total degree of / and p kt i is a univariate polynomial of 
degree < deg f — k. This representation is often called the Gauss represen- 
tation. A similar formula is valid for a much larger class of functions. Let us 
recall that a function / : G — > C defined on an open set G in M. d is called 
polyharmonic of order N if / is 2N times continuously differentiable and 

A N u(x)=0 (1.9) 

for all x e G where A = + ... + is the Laplace operator and A N the 
N-th. iterate of A. The theorem of Almansi states that for a polyharmonic 
function / of order N defined on the ball Bjj = {xeR d :|x|<i?} there 
exist univariate polynomials p k ,e (r) of degree < N — 1 such that 

OO Qfc 

/(*) = £X>M(M a ) y M(*) ( L1 °) 
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where convergence of the sum is uniform on compact subsets of Br, see e.g. 

In this paper an important role is played by the representation of poly- 
harmonic functions on the annular region A a ^ developed by Vekua and 
Sobolev (see [54], [47]): It is especially simple if the dimension d of the eu- 
clidean space K d is odd. Then for any polyharmonic function / of order N 
defined on the annular region A a ^ there exist univariate u k ,i (•) , u k ,i (•) poly- 
nomials of degree < N — 1 such that 

oo a k 

«W=EE («W (M") rk + «M (|x| 2 ) r - d - k + 2 ) Y klt (9) (1.11) 

where convergence of the sum is uniform on compact subsets of A a ^. An 
equivalent way to express this is to write 

(N-l N-l \ 

Y^W k+21 + £<W "' Y k , e (9) . 

3=0 3=0 / 

In particular, one obtains a representation of the polynomials if all d k ^-j are 
zero. 

1.2. Complexification of the annulus and the annulus of the Klein-Dirac 
quadric 

We want to study analytical extensions of functions / defined on the annular 
region using the Laplace-Fourier series. Instead of working with functions / 
which are a priori analytically extendible to a fixed domain U in the complex 
space C d we shall require only that we can extend the function x = r9 i — > 
f (r9) to an analytic function z9 i — > f (z9) , so we only complexify the radial 
variable r to a complex variable z. Then one should expect from (|1.7I) that 
the Laplace-Fourier coefficient fk,£ (r) extends to an analytic function of one 
variable. 

To make things more precise let us introduce the basis functions 

Rj W := ( r -d-k+2+2U-N) for j = JV.JV + 1, 27V - 1. (L12) 

By representation p. lip we sec that the functions 

{Rj (r) Y k/ (9) , For k > 0,1 = 1, 2, a k , j = 0, 1, 2N - 1} 

form a basis for the functions polyharmonic of order N in the annulus A a ^ 
For every fixed k we introduce the following convenient notation for the 
set of indices to be used further: 

lZ k := {j E Z : j = k + 2m, or j = —d — k + 2 + 2m for an integer m > 0} . 

(1.13) 

As said above, the main feature of our approach is to work with a com- 
plexification of the functions / defined in the annulus A a ^ of the form / (z9) . 



lr The functions Rj for j = 0, 1, .., N — 1 are called polyharmonic Kelvin transforms of the 
functions ijjy+j, cf. [3], p. 6. 
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For polyharmonic functions / representation (jl.lip provides a direct descrip- 
tion of this complexification by considering the complexified basis functions 
Rj (z); we consider the linear space generated by the functions: 

V k>d :=span{{z 2 ^}; =0 , K^ 2 ^'}^} (1.14) 

Now in these new notations, a generic element in the space Vh,d has the 
form Uk,e (z 2 ) z k + u ki £ (z 2 ) z~ d ~ k+2 where u kj ( (•) and u k j (•) are polyno- 
mials. Obviously, V k , d 'c H (A 0>6 ) , V k ,d C H (A , 6 ) , and V k , d C H (C\ {0}) 
where "H (f2) denotes the space of all analytic functions in a domain fi c C. 

Thus, we have come naturally to the next objective of the present re- 
search, the Function theory in the space C x S™ -1 . The main interest in the 
present paper is devoted to the complexified annulus A a ,b given by 

A a ,b ■■= {{z,0) ■ a < \z\ < b, 6e S^ 1 } = A a , 6 x S^ 1 ; (1.15) 

here A a! b is the annulus in C, 

A a ,6 := {z e C : a < \z\ < b} . 

Our main novelty will be a multivariate generalization of H 2 (A Qj b) 
called (polyharmonic) Hardy space on the annulus A a ,b to be introduced 
in Definition 14.51 We will denote this space by H 2 (A a ,b) since it depends 
on a parameter L > 1. The name polyharmonic comes from the fact that 
H\ (A a .b) is obtained as a limit of the complexifications of the polyharmonic 
functions in the annulus (jl.lip . namely, of the space: 

V := span {z fe + 2 ^ M (6) , z - d - k + 2 +^Y Kl (9) : k, j > 0,1 = 1, 2, a k ) , 

(1.16) 

where z £ C and 9 E E> d 1 . Clearly, the elements of V are finite sums of the 
type 

u (z, 9) = J2 («M (z 2 ) zk + (z 2 ) z- d - k + 2 ) Y k j (9) , (1.17) 

k,t 

where u k ^ (•) and u kt i (•) are algebraic polynomials. In view of the Poly- 
harmonic Paradigm announced in [3D] , the space H 2 (A a ,b) generalizes the 
classical Hardy spaces which are obtained as limits of (Laurent) polynomials, 
where the degree of a polynomial is replaced by a degree of polyharmonicity. 

Let us summarize briefly our objectives: The main approach to the 
polyharmonic Hardy space H 2 (A a ,b) is to generate it by the complexified 
finite sums (|1.17p . by introducing an appropriate inner product (and norm). 
This Hardy space will be a Hilbert space and we will provide a Cauchy 
type kernel, which is the analog and a generalization to the Hua-Aronszajn 
kernel in the ball (cf. 3 , p. 125, Corollary 1.1). Let us note that the last 
is a multidimensional generalization of the classical Cauchy kernel - from 
Complex Analysis (various about Cauchy kernels see in [H], |28j). 
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Remark 1.1. Our results are closely related to the Function theory on the 
Klein quadric (called sometimes Klein-Dirac quadric) given by: 

KDQ := {z9 : z £ C, G S^ 1 } =Cx § d - 1 /Z 2 , (1.18) 

which may be considered as a subset in C d . The factor Z2 means that in the 
set C x S d_1 we identify the elements (z, 9) rj (— z, —&) and the elements zQ = 
(z9i,z0 2 , zOd) of KDQ are points in C d . The space KDQ is very important 
for us since a genuine counterpart of the present Hardy spaces exists only on 
the ball B x § n - 1 /Z 2 of KDQ (and not on the ball © x S"" 1 ) which is due to 
the fact that the polyharmonic functions in the ball in IR d have representation 
of the type (jl.lip with terms z k+2j Y ky i (9) = z 2l Y k ^ (z9) depending on z9. 
In a forthcoming research [26], [24] we study Hardy spaces on the ball in 
KDQ . For the Klein quadric the ball D x § n_1 /Z2 is the compactified four- 
dimensional Minkowski space-time S 1 x § 3 /Z2. It is hoped that the present 
research provides a useful framework for CFT by approximating the fields by 
elements of the Hardy space in the annulus A a ,b, see Remark 14.91 below. 

The Klein-Dirac quadric KDQ defined in (|1.18[) for arbitrary dimension 
d, has been originally introduced in a special case by Felix Klein in his Erlan- 
gen program in 1870, where he put forth his correspondence between the lines 
in complex projective 3-space and a general quadric in projective 5-space. The 
physical relevance of the quadric and the relation to the conformal motions of 
compactified Minkowski space-time had been exploited by Paul Dirac in 1936 
[12] . The Klein-Dirac quadric plays an important role also in Twistor theory 
where it is related to the complexified compactified Minkowski space 38] . 
Apparently, the term "Klein-Dirac quadric" for arbitrary dimension d, has 
been coined by the theoretical physicist I. Todorov, cf. e.g. [35], [35]. In these 
references important aspects of the Function theory on the ball in KDQ have 
been considered in the context of Conformal Quantum Field Theory (CFT). 
In the context of CFT Laurent expansions appear in a natural way as the 
field functions in the higher dimensional conformal vertex algebras (using a 
complex variable parametrization of compactified Minkowski space); see in 
particular formula (4.43) in [37] . as well as the references [33], [36] . 

The paper is organized as follows: in Section [5] we recall background 
material about the Hardy space H 2 (A a ^) . In Section [3] we define the "com- 
ponent" subspaces H 2 ' k (A a- &) of H 2 (A a ^) which will be necessary for the 
multivariate case. In Section [3] we introduce the polyharmonic Hardy space 
H\ (A a ,b) on the complexified annulus A a ,b- We prove that it is a Hilbert 
space, a maximum principle, and infinite-differentiability of the functions in 
H\ (A a ,b) ■ In Section [5] we construct a Cauchy type kernel for the compo- 
nent spaces H 2 ' k (A a ^) and in Section El a Cauchy type kernel for the space 
H 2 (A a ,b) ■ In Section[7]we prove the main result of the paper, about the error 
estimate of the polyharmonic Gaufi-Jacobi cubature formulas introduced in 
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2. The inner product and norm in the annulus in C 

Let us recall the classical Hardy space H 2 (A 0l b) on the annulus A a ^ and its 
basic properties, cf. [43], p. 4; see also [5]. 

Let us remind the usual mean value of a function / defined in A a .t, as 



Definition 2.1. The Hardy space H 2 (A aj b) is defined as the space of all ana- 
lytic functions / in A a> b by the condition 



Remark 2.2. Recall that by a theorem of Hardy and Riesz the function 
log M2 (/; r) is a convex function of log r, hence its maximum is attained 
either for r — > a+ or for r — ¥ b— ; cf. |43j . p. 4. 

The following Theorem contains the basic properties of the space 
H 2 (A a , b ) , see 03], p. 7-9. 

Theorem 2.3. Every f € if 2 (A^f,) is decomposed in f = g + h where g £ 
H 2 (D b ) and h e H 2 (C\ 57) , w/iere D r denotes the disc {z G C : |z| < r} . 
(Hence, the functions g and h provide the limiting properties to f ) 

Let f € -ff 2 (A 0l b) and /* oe its limiting function on dA ay i, which we 
denote by f* on the circle a x S 1 C C, and by / 6 * on £/ie circle b x S 1 C C. 

(%) If the Laurent expansion 0/ / is given 6j/ 




(2.1) 



sup M 2 (/; r) < 00 



(2.2) 



a<r<6 



00 




(2.3) 



then faifb € L2 (S 1 ) and /iave Fourier series respectively 





(2.4) 



Vace versa, to every pair of functions fa,fb ^2 (S 1 ) having Fourier series 
\2.J$ there corresponds uniquely an f G i/ 2 (A a) b) given 01/ if£.ff|J . 
fii) The following Cauchy formula holds: 



1 f 27r be i{ f 1 /" 27r np iL P 



(Hi) The following Maximum principle holds: 




for all a < \z\ < b. 



8 



Ognyan Kounchev and Hermann Render 



Remark 2.4. The space H 2 (A a fi) is a Hilbert space with the inner product 

!l/llff 2 (A a ,t,) = ll/olli 2 + ll/tllia 

1 r 2lx 1 r 2lx 

= lim — / f (re lv ) g (re l v)d(p + lim — / f (re ltp ) g (re ll ?)dtp 

r^a2TTj r >b 2n J Q 

cf. [33]. Let / £ H 2 (A aj b) and have a Laurent expansion 

oo 

f(z)= J2 ^. 

3=-oo 

Then 

OO OO 

\\f\\W a ,) = E 1*1 VllL (< *. lt) = ^ E i«*i a (a«+^). 

j=—oo j=—oo 

Remark 2.5. In [7] S. Bergman computes the kernel for the annulus, cf. also 
[2J and [32]. 

Remark 2.6. The following follows directly from above: For every function 
f E H (A a ,b) we have the inequalities 

ll/ll ff2(Aoi6) < 2 sup M 2 (/;r)< 2 ||/|| ff2(Aa 6) . 



3. The component one- dimensional spaces H 2 ' k 

For us important role will be played by the so-called "component Hardy 
spaces" H 2 ' k (A a ,b) which are relatively simple subspaces of the Hardy space 
H 2 (A a> b) , cf. some classical references on Hardy spaces as [T7], [4"U] . 

Definition 3.1. The closure of the linear space of analytic functions Vk,d de- 
fined in (|1.14p . in the Hardy space H 2 (A a .b) will be denoted by 

H 2,k = H 2,k (A a ,b) (3.1) 

and will be called component Hardy spaces in the annulus A Q) b. 

The properties of the space H 2 ' k (A 0j b) easily follow from those of the 
usual Hardy spaces H 2 {A a ^) and H 2 (B&) . 

Proposition 3.2. Let f E H 2 ' k (A a . b ) . 

1. Then z d+k ~ 2 f (z) € H 2 (B b ) . 

2. There exist two uniquely determined functions f\, fi G H 2 (Bb) such 

that 

f(z) = z k f 1 (z 2 )+z- d - k + 2 f 2 (z 2 ). (3.2) 

3. Let f* (having two parts, f* and f b * ) be the limiting function for f 
on dA at b. Then the following expansions hold, 

f a (oe*) = c J ajeijV > fb (be*) = E ( 3 ' 3 ) 
jeiz jell 
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where the set 1Z is defined in M.13\) . Vice versa, for every two functions 
f* (ae lv ) and f* (be iv ) with the above expansions there exists a function 
f 6 H 2 ' k (A a having Laurent expansion V2.3\) for which f* is a limiting 
one. 

4- There exist constants A,B>0 such that for every f £ H 2,k (A 0j b) 

holds 

A \\f\\w(A a , b ) ^ ll/lllfP(Bi,) + ||/2||_H-2(D b ) ^ B II/IIh 2 (A q>(i ) ; 

cf. 03], p. 5. 

4. Polyharmonic Hardy space H\ (A a ,b) on the complexified 
annulus 

4.1. A motivation 

Now we may define the polyharmonic Hardy space i?£ {A a , b) on the com- 
plexified annulus A a .b of R d . The most important property of this space will 
be a Maximum principle, and a related Cauchy type kernel. 

Let us explain the problem of finding a proper norm to define the space 
H\ (A a ,b) ■ The elements of the space V, which is generating if£, are functions 
/ (z,9) having finite expansion 

/ (z, 9) = £/ M [z) Y k j (9) = ]T (f k j (z 2 ) z k + J M (z 2 ) z - d - k+ A Y kJ (9) . 

k,t k,t ^ ' 

(4-1) 

Since all components of this sum are mutually orthogonal in 9 E S d 1 , it 
would be natural to consider as a candidate the norm 

11/11* : = X! \\fk/\\ H 2,h ■ 
k,e 

However it is easy to see that such a norm does not guarantee that a Cauchy 
sequence of elements /at (z9) , N > 1, will converge uniformly on compacts 
of A a ,b- Indeed, let us consider the following sequence, 

u N (z,0) = z -d-k N+ 2 + 2N Yk N A9) = Y kN A9) (4 2) 

\Jk N \Jk N 

with fcjy = — d + 2 + 2N, for N > d. Obviously, ujy £ V and we have 
||« w ||, = for TV 

However, as we see below this sequence is even unbounded (for d > 3 ). 

Proposition 4.1. For d > 3, i/ie sequence of functions un defined in \4-2^ is 
unbounded. 

Proof. First of all, recall that for every spherical harmonic Y k ,i (0) of degree 
k > (sec (11.41) ) we have the estimate 

\Yk,e (0)\ < Cki- 1 for all 9 e S d_1 . (4.3) 
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(cf. [IS]). Let us consider d — 3. We have the following representation of the 
spherical harmonic which corresponds to the Legendre polynomial P k (•)> (cf. 
[55] . Appendix 2, part II ): 



2k 4- 1 

Y k i (6>) = \ —r—Pk (costf) for < < tt, < ip < 2tt; 

V 47T 

here we have 9 = (cos <p sin sin <p sin cost?) € S 2 , and the usual Legendre 
polynomials Pk (•). Since Pk (1) = 1, it follows that the sequence ujv ( z ,9) is 
unbounded at # = 0. □ 



Remark 4.2. An additional observation about the unwished behavior of the 
sequence itjv {z, 9), which is easy to see for d = 2 (where Y k) \ = cosfctp and 
^fc,2 = sin kip) is that derivatives in 9 of ujv (2, 0) contain non-negative powers 
of'fc. 

For that reason a proper norm has to "tame" the behavior of such se- 
quences including their derivatives: we will introduce an inner product and 
norm which contain weight factors for every component k, which "suppress" 
such unwished behavior and this will guarantee convergence of Cauchy se- 
quences in H 2 (A a ,b) on compact subsets of A a ,b- 

Remark 4.3. In this context let us remind that a function / (9) on S d_1 

00 

is real analytic on if its Laplace-Fourier expansion / (9) = V^^fc (9) 

k=0 

a-k 

(where Y k (9) = }Jk,tY kii (9) ) satisfies 

i=\ 

\\Y k (d)\\<Ce- r > k forfc>l 
for some constants C > and i] > 0, cf. [47] . 

4.2. The weighted inner product and norm 

Let us introduce the weight factor to be an arbitrary number 

L > 1. 

Definition 4.4. Let / := {fk,i) k t an d 9 '■= {9ki) k t be two sequences of 
functions with f k> t e H 2 ' k (A a ^) and g k ,£ G H 2 ' h (A 0j b) for all k = 0,1, ...,and 
I = l,..,a k - We define the following inner product 

</,<?) ff = ■■=Y,(f^9k,i) HHAab) L 2k , (4.4) 

k,e 

where (,)#2( A b \ is the inner product in H 2 (A a> i,) . 
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It is easy to see that (,) H 2^ A b \ satisfies all necessary properties of inner 
product on a Hilbert space, and the norm is given by 



ii/ii V= ii/ii^(^) : =Eii^ii^(A. fi ) ia * ( 4 - 5 ) 




r e ^ r ^ ^ i 2 ^ + L2k r ^ i 2 ^ 



Let us consider the function 



with /fc^ 6 i/ 2,fc (A a ,b) and ||/||#2 < oo. By applying Proposition 12.61 we see 
that for every r G (a, 6) holds 



which shows that / is a function in Li (A a .t) ■ 
Now we may give the following definition. 

Definition 4.5. The polyharmonic Hardy space H\ (A a ,b) consists of all func- 
tions 



with f k .e e H 2 k (A a ,b) and satisfying ||/|| ff 2 < oo. 

Let us remark that this definition is completely analogous to the classical 
case where one takes the Laurent series of the type J^^_ 00 a,jZ J which are 
convergent on compacts in the annulus A a .fc. 

It is easy to see that there are equivalent representations of the functions 
in the space H\ (A a ,b)- 

Proposition 4.6. Every element f € H 2 (A a ,b) has the following representa- 
tion as an infinite series (in L 2 — sense): 




2 



f(z,6) = J2h,t (z)Y ktt (9) 




(4.6) 
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where for every index (k, £) the function f k ,e (z) € H 2,k (A 0j b) ; the set IZk 
was introduced in \1.13\) . The norm of f is given by 




h,:/(<r ! ■ lr J ) • 
z„ J 

The proof is straightforward. 

We prove the following fundamental result which contains as its first 
item that the space H 2 (A a ,b) is a Hilbert one. 

Theorem 4.7. 1. The space H 2 ^ (A a ,b) "is complete. 

2. The following Maximum principle is true: for every compact K C 
A a ,b there exists a constant C > such that 

< ■ tz | | | , ttII/II^ca,,) f° r i^)€K. 

mm (1 — \z\ /o, \z\ /a — 1) ^^^^j 

S.Iffe H 2 (A,b) then f e C°° . 

Proof For the proof of item 1), let {f n } n>1 be a Cauchy sequence in (A a ,b) ■ 

Then it is clear that for every index (fc, £) the sequence of functions I fj} e (z) > 

I ' J n>l 

is a Cauchy sequence in H 2 (A aj &) (even in H 2 ' k (A a ^) ), and as such has a 
limit g kA (z) € ff 2 ' fc (A a . b ) . 

We have to show that the function defined by 

g(z,e) = J29k,t {z)Y k<t (9) 
k.e 

belongs to H\ (A a ,b) and is the limit of {f n } n >i there. Let us fix some e > 0. 
Let us take some k\ > 1. By the Cauchy sequence, we find rii > 1 such that 
for all n, m > n\ we have 

fcl a k ^ 

ZJIJ II^M ~ •^llff 2 (A a , i) ) L2fc - II/" ~ f m WHl(A a , b ) < 

Taking the limit for m — > oo we obtain 

\\fk,i - 9k,e\\ H *(A atb ) L2k < £2 - 

k=oe=i 

Then taking the limit k\ — > oo ends the proof. 

2). Let us consider an annulus A a > y with a < a' and b' < b which 
satisfies K C A a '.b'- For a function / given by (|4.1j) . by the triangle inequality 
we obtain for all (z,9) £ K the following inequality: 



\f(z,e)\<J2\fkA*)\\Yk,i (0)1 . 



fc.f 
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By the maximum principle for the annulus, Theorem 12.31 and the estimate 
for the spherical harmonics (|4.3p we obtain 



^ * min(lH^N/ fl -l) g M 

Further, by applying the Cauchy-Bunyakovsky-Schwarz inequality, for every 
(z, 0) £ A a ',b' we obtain 

min (1 - \z\ /b, \z\ /a - 1) 



2^ ^ d 

- minfl-|z|/6Jz|/a-l)^ Lfe "-M^OW) 

1/2 / \ 1/2 

2C 

< 



min(l-|z|/6,|z|/ 



2 A- 



Since the constant in (|1.4j) is estimated by = O (fc d 2 ) (cf. [SO], chapter 
4.2) it follows that for all (z,9) G holds 

I"**)" ^ min(l-|^N/ a -D " / "^^ X (p d - lkd - 2L - 2k ) 1 



< 

mm ( 1 — 



\k=0 

C * II f \\ 

(l-\z\/b,\z\/a-l) llJllH K A °^ ■ 



This ends the proof. 

3. As above we consider A a \b'- Let / £ V have the expansion 

f(z,e) = J2fk,e {z)Yk,i (0)- 

k,t 

For (z,6) £ A a ',b' we obtain the estimate 

\df( z ,e) /dz\<^2\fi, oolm,* wi- 

By the usual Maximum principle in the annulus A at b (which is equivalent to 
continuity of point evaluation operator in H 2 (A 0j b) , cf. section 17.8 in [4"U] ) 
we obtain 

\df(z,e)/dz\<cJ2\\fkA H ^A a , b ) ki ~ 1 fOTa11 (z,0)eA a%v , 

k,t 

and the last is convergent as we have seen above. In a similar way, for some 
derivative Dg in 8, where a is a multiindex, we apply the following inequality 

\D%Y k ,e (0)| < Cifcl Q l +i ^ for k > 1, 9e § d ~\ 
cf. gS], p. 120. This ends the proof. 

□ 

Remark 4.8. Formula (|4.6|) is a generalization of the Laurent expansion. 
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Remark 4.9. Let us indicate a relation between the Hardy space H\ (A a ,b) 
and analogous space on the ball D x § d_1 /Z2 of the Klein-Dirac quadric KDQ 
given in (|1.18l) . Let the function / e H\ (A a ,b) be polyharmonic of order N 
in the sense that in representation (|4.1[) 

/ (z, 9)=J2 (fk,i (^ 2 ) * k + hi (z 2 ) z- d - k+2 ) Y k4 (9) 

k,l ^ ' 

the functions fk,e, j ' \. i are polynomials of degree N — 1. Then we have the 
representation 

f(z,8) = f 1 (z,6) + z- d+2 f 2 (z,9), 

where the functions 

h{z,6)=YJk,l (z 2 )YkJ~ z 
k,e ^ z 



a obviously defined on the KDQ \ {0} . 



5. The Cauchy type kernel in the component spaces H 2 ' k (A aj &) 

By the Cauchy formula in the annulus A a ^b, for every function / which is 
analytic in a neighborhood of A^f, we have: 

OO j 



Jr b t — z 2ni J Ta t — z 



2lT ./n I T J 



I) 



dp, 

T =ae' , P 



where we assume that the circular contours T a = a x S 1 and = b x S 1 have 
the "positive" (anti-clockwise) orientation and are parametrized by t = ae %v 

and r = be lip ; we used the notation g (r) — g (be lv ) . Hence, we select 



out of the above series only those terms which belong to the space (|1.14[) and 
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we define the following kernels: 



k+2j 



m=-d-k+2+2j>0 



for \t\ = b 
for |t| = b 



Kl(z,r):= J2 (-) for|r|=o. 



(5.1) 
(5.2) 
(5.3) 



m=-d-/c+2+2j<0 

It follows that for every / € Vfc,d holds 



f(z) 



1 

2W 



(K 1 (z,r)+K 2 {z,r))f(T) 



T=be i f 



1 f 2 " 

tJ ^(*,r)/(r) 



dtp 



Now for every integer fc > we will define the Cauchy type kernel for the 
annulus K k by setting: 



K k (z, t) := K\ (z, t) + K 2 k (z, r) for |r| = & 
K k (z, r) := #| (z, r) for |r| = a. 



(5.4) 
(5.5) 



We have the following Proposition which shows that the function K k 
is the Cauchy type kernel for the space H 2 ' (A a! b) endowed with the inner 
product in Remark [ 



Proposition 5.1. Let us denote by f* the limiting value of the function f G 
H 2 ' k (A a ^) on dA a ^. Then for all z 6 A a .f, the following Cauchy type formula 
holds 



1 



2^ Jo 



f(z) = — (Kl(z,T)+Kl(z,T))r(r) 



T=be i! -P 



dip 



i- Kl{z,r)r(r) 
2lT Jo 



d(p 



1 



K k (z,r) f* (r)-rfr 

dA a . b T 



= (Kk(z,-),f(-)) H s(A a , b ) ■ 

The proof follows directly from the Cauchy formula for the Hardy space 
H 2 (D c ) (see also Theorem 1 on p. 9 in g3]). 

A nice feature of the kernel K k is that we have concise expressions: 
Obviously, we have 



ft x 2 
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For k odd we have the expressions 

00 9 1 2 

*2 = E = izk* = ^ for |t| =6 



p=0 

ri + fc-2 



-(f)' 



- /T x2p / T d+fe-2\ 



For fc even we have the expressions 



^(^-) = E(7) 2P+1 = 7CT fOT H = 6 



p=0 

d+fc-2 



K A k (z,r)= > (-] = — |1- „ ) x Ki(z,T) for in = a 



3 I*) ' ) - /_ 

=.1 

The following estimates are important. 



Theorem 5.2. For every e > wif/i £ < (6 — a) /3 and /or z and r satisfying 
a + e<\z\<b — s and |r| = a, or |r| = 6, and /or all k > holds 

\K k (z,r)\<C £ , 

where the constant C e > is independent of k. 

The proof follows directly from the above expressions for the kernels 
#£,.7 = 1,2,3. 



6. Cauchy type kernel for the complexified annulus A a ,b 

Finally, we are able to construct the Cauchy kernel for the complexified an- 
nulus A a .b- 

First, we prove the following: 

Proposition 6.1. Let f £ i?£ (A a ,b) ■ Then in L2 sense the following limits 
hold true: 

lim / (re*, 6) = ( ae ^) F M (*) 

r—ta z — * ' 
k,£ 

lim / (re*0) = (be*) Y k/ (6) . 

k,t 

The proof mimics the classical proof in [JD], Theorem 17.10. 
Proof. First of all, we see that the function 

g(e*,B)=Yjh (^^C) 
k,i 
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2tt 



J^fli (afi V ) Y kA (9) 



k.t 



d9d<p <J2\\tk,e 



L 2 (dA atb ) 



k,l k,i 

< Cl ll/llff£(,A 0ii) ) > 

where C, Ci > 0. Similar estimate holds for b. 
Let us put 

^ ; = / r ( ae ^) y M - e/m (^^) r M 

We have to prove I r — > for r — > a. We have obviously 

,-2n 



dOdip. 



/>Z7T 

Ir < J2 / |/M (° e ^) " /M (^) I 

We fix some e > 0. There exists a fci > 1 such that for every k > k\ 



holds 



This implies 



EE H.MIff2 (Aoi!)) I' 2fe < |- 



k=k 1 i=l 



00 a fc <.2ir 

E E / l/M ( ae ^) ^ /M (re^)|% < C| 

for some constant C > 0. On the other hand, we find 8 > such that for 
Ir — a| < 5 holds 



fel— 1 ifc r 2ir 



E E / I/m ( aelip ) - /m (™ ,v ) r ^ < c-. 

Hence, we obtain / r < Ce which proves that I r — > for r — > a. 



□ 



Recall that for every / e H\ (A a ,b) all components fkj S _ff 2,fc (A 0j b) 
and they have limiting functions f^ e £L 2 (dA a: t) ■ Hence, by Proposition ^. II 
we obtain 

/ (z, 9) = (Kk (z, ■) , fZ ti (-)} HHK b) Y k , e (9) for all (z, 9) e A a , b . 

k,e 

Definition 6.2. We define the Cauchy type kernel K for the space H\ (A a ,b) 
by putting 

oo a k 

K (z, 9- r, 9') = E Lk K * ( z ' T ) *M (*) r M • (6- 1 ) 

fc=0 £=1 
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The function K (z, 9; r, 9') converges on compact subsets of A a ,b'- Indeed, 
we use the estimates |Yfc,^(0)| < C\k~^~ , and a k < C2k d ~ 2 , cf. [45]. By 
Theorem 15.21 we obtain the estimate 

oo id—2t.d—2 

\K {z, 9; r, 6')\ < C 3 , £ Lk < oo 

fe=0 

for all a + e < \z\ < b — e and |r| = a or \t\ = b. 

Finally, we generalize the classical one-dimensional Cauchy formula, and 
justify the name "Cauchy type kernel" given to K (z, 9; r, 9') in Defmition. 16.21 

Theorem 6.3. The kernel K (£, 9; z, 9') is a Cauchy type kernel for the space 
H\ (A a ,b) defined by the inner product h4-5\ ), i.e. for every f g H\ {A a ,b) 
holds 

(K (z, 9; z', 9') , f* (z\ 6')) Hl = f (z, 9) , for all (z, 9) E A a , b , (6.2) 
where f* is the boundary limit at dA a ,b of the function f. 
The proof follows by a direct application of Theorem [ 



7. Error estimate of the Polyharmonic Gauss-Jacobi Cubature 
formula in the amiulus 

The topic of estimation of quadrature formulas for analytic functions is a 
widely studied one. Beyond the classical monographs [3T], [TT], we provide 
further and more recent publications, as jS], [H], [H], [TH], [27], [S3]- No 
references may be found though for the the Numerical Integration in the 
multivariate case, which is often called cubature formulas, even in the funda- 
mental monographs as [5T|, @H]; see also the recent survey [5]. 

In order to explain our approach, first of all, we will recall the polyhar- 
monic Gauss-Jacobi cubature formula presented in [21] . [23] for the annulus 
A a j, C R d . Let us fix some integers N > 1 and k > 0. We consider the 
following 2N— dimensional subspace of Vk,d in (|1 .14[t . 

V k , d , N := span^r^to 1 - {r^^ }^} • 

We remark that this is a Chebyshev system of order 2N, cf. [50] . 

Let the pseudo-positive (signed) measure d[i be given in the annulus 
A a ,b, see (|1.5p . For all indices (k,£) the component measures are defined by 

dnwjt (r) := / Y kl (9) d^ (r9) > for all r e [a, b] ; (7.1) 

here the integral is symbolical with respect to the variables 9. Rigorously, the 
component measure dfik,e (r) is defined for the functions g (r) on the interval 
[a, 6] by means of the equality 

/ 9 (r) dnk,t (r) := / g (r) Y kJ (9) dfi (x) ; 

cf. E5, US]- 



Hardy Spaces and Cubature Formulas 



19 



Let us fix (k,£) . Recall from [5T], [23], that we consider truncated 
Moment problem and related Gauss- Jacobi type Quadrature formula corre- 
sponding to the operator A 2N . Thus, there exist points tk,e-,j, j — 1, 2, 27V, 
belonging to the interval [a, 6] , and non-negative numbers {^k,£-.j}j =1 , such 
that the following quadrature formula holds: 

h rb 

Q (t) d/* M (t) = Q (t) di% t (t) (7.2) 



= /M,e-,jQ (tk,£;j) for every Q £ V k ,d,2N, 

3 = 1 

cf. [30], Theorem 4.1, chapter 4. We have put 

N 

i=i 

First of all, we will find the error for the above Quadrature formulas 
(|7.2[) following the classical scheme outlined in [TT], p. 231 — 235 (see also 
chapter 12 in |31j). 

Let Hk,t,2N [/] (t) be the unique element of Vk,d,N which interpolates 
the function / at the points tk.e-j, i.e. 

Hk A2N [/] = / (*Ma) for j = l,2,...,2iV. (7.4) 

On the other hand, for every function / € H 2,k (A. a ,b) by Proposition 15 . 1 1 we 
have the Cauchy type formula 



/(*)=«^/ K k (z,r)r(r)-dT, 

2m J dknh T 



where /* is the limiting value of / on the boundary dA a ,b- Hence, we obtain 
H k ,t.2N [/] (z) = / H kA2N [K k (•, r)] (z) f* (r) -dr; 

here for every fixed r, by H k j-2N [K k (-,t)] (z) we denote the element in 
Vk,d,N which interpolates the kernel function K k (z, r) at z = t k j.^ for j = 
1,2, ...,2N. This gives us the remainder formula of Cauchy type for the in- 
terpolation ([73 . 



/ (z) = H k . e ,2N [/] (z) (7.5) 
+ I {Kk (z, t) - H k j, 2N [K k (; t)} (z)} f* (r) -dr. 



Now we are able to estimate the error functional: 

rb n b 

,G 



E k ,e [f] := / / (t) dnw (f) - / / (t) dft. (t) . (7.6) 
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Indeed, by (|7.5I) we have 



/ (z) dnk,l (z) 



1 



H k ,£-2N [/] {z) dfik,e (z) 



{K k (z, t) - H k ,t;2N [K k (•, r)] (z)} d/i M (z) /* (r) -dr 



On the other hand, by (|7.2I) and (|7.4I) we obtain 



H k ,i-2N [/] (-2) ^Mfc^ ( z ) 



•ffk,/;2iV [/] (z) d/iw (z) 



f(z)dfi? e (z) 



which implies 

E k4 [/] = / / (t) d/x M (t) - I f (t) dfi. (t) 



(7.7) 



= tt~ / / {^fe t) - H kl , 2N [K k (•, r)] (z)} d MM C*) /* ( T ) 

Zm JdA a ,bJa T 

and we see that in the right-hand side we do not have derivatives of the 
function / but only its values on the two circles <9A a .t>. Now we may prove 
the following Lemma. 

Lemma 7.1. Let the knots {tfe,^ } 2 =1 of the quadrature formula J7.ff| ) lie in the 

interval [a, 6] C (a',b') . Then for every function f in the component Hardy 
space H 2,k (A^b') , the error of the quadrature formula {7.2ty satisfies 

\E k Af]\-- 



/(i)d/i M (t)- / f{t)d^At) 



<Ck\\f\\ H , iAa , b ,y (7. 



Proof. Since {tk,i;j} 2 = i € [«, 6] C (a! ,b') it follows that the interpolation 
operator H k ,i;2N [Kk (•)!")] depends continuously on r g dh^p. Hence, the 
function 

/ {K k (z, t) - H kA 2N [Kk (•, t)] (z)} d^M («) 

^ a 

is also continuous in t S <9A a <.f/ and by Theorem [521 its modulus is bounded, 
namely 

rb 



/ {iffc (z, t) - H kA 2N [K k (•, r)] (z)} dfi k ,i (z) 

./ a 



<c fc 



for t 6 



By application of Cauchy-Bunyakovksi-Schwarz inequality, from (|7.7p . we 
obtain for some constant C > the inequality 

< CC k 11/11^2^,^,) 



f(t)d^(t)- I f(t)d^ e (t) 
which ends the proof. 



□ 
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Finally, we may proceed to the polyharmonic Gauss-Jacobi Cubature 
formula for the pseudo-positive measure d/i. We will consider this in a simpler 
setting, for functions / having a finite Laplace-Fourier expansion: 



/(*)=/ (r6 



J2 fk,t(r)Y k<t (9) 



(7.9) 



k,l, k<k 



The polyharmonic Gauss-Jacobi measure d\i (x) is defined as the pseudo- 
positive measure having component measures (|7.1|) equal to f in (|7.3p . 
Written symbolically, we have the expansion 



dfi G (x) = 5> M (0) dfj%i (r) . 



(7.10) 



kj. 



The corresponding Cubature formula is simple to formulate if the function / 
defined in A a & has a finite expansion as in (|7.9|) : 



£ [ b fk.Ar)d^ e (r 



(7.11) 



It is direct to see from the quadrature formulas (|7.2[) . that for a function u (x) 
which is a finite sum of the type (fTTTj) with A 2N u (x) = in A a<b , holds 



u (x) d[i (x) 



u (x) d/i (x) , 



cf. [21) . |23) . This justifies the name polyharmonic cubature formula. Hence, 
we define the Cubature formula by putting 



CnU) 



f{x)dy G (x). 



(7.12) 



We will prove the estimate of the cubature formula in the following. 

Theorem 7.2. Let d/i be a pseudo-positive measure, and Cn (/) be the poly- 
harmonic Gauss-Jacobi Cubature formula defined by in the annulus 
A a ,b- Let a 1 < a and b < b' , and we assume that the function f € U\ (A a > ,b') 
has a finite representation |7.ff[ ) for some kg > 0. 
Then the error of the Cubature formula J 7. 1 l\j is estimated by 



\E[f]\ = 



f (x) dfi (x) - C N (f) 



< 



CI 



> E rjk x Wf\\ 2 Hl(A a ,y)- 

\ k,t, k<k 



Proof. Due to the finite expansion, we obtain the equality 
f (x) d[i (x) = 

k,t, k<k " 



^2 / fk,e( r ) d li-k,i(r) 
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On the other hand, by the Cubature formula (|7.11[) we have 

f f(x)d^ G (x) = fhAr)d^(r) 
jA ". b k,l, fc<fc 

Hence, we obtain the error functional 

E[f]= f f(x)dn(x)- f f (x) d/i G (x) 

JA a , b JA a _ b 

cb 

fk,i (r) dfi G j, (r) 

k,£, k<k a \ Ja 

and by Lemma 17. II we obtain 



^2 / fk,t M d^k.e (r) - / 

L k<k„ \ Ja Ja 



\E[f}\ = 



f{x)d l x{x)- f(x)d^(x) 

A a . b JAa,.b 



k,i, k<k a ' k,l, k<k 



< 



\ 



C 2 

Y J2k x WfWm 



k,l, k<k 

This ends the proof. 



□ 



For practical purposes it is reasonable to work with functions having 
finite expansion (|7.9[) . However, from Theorem l7.2l we see that the asymptotic 
behavior of the constants Ck defined in (|7.8JI is important for the estimate of 
the error functional E [/] for functions / which do not have a finite expansion 
(|7.9|) . This is a subtle task of Interpolation theory which is beyond the scope 
of the present paper. 

Another important remark is that the cubature formula (|7.12|) may 
be easily extended, although in a non-unique manner, to signed measures 
having bounded variation. Indeed, we may represent the measure as 
/i 2 , where (j, , ji 2 have bounded variation. We can 



a difference, 



Ml. 



4,. 



do this representation componentwise (non-uniquely) by \Xk 
with measures d\i\ d^ 2 ^ > 0. Then we find the poly harmonic Gaufi-Jacobi 

cubature measures and ^ and put C N (/) = / fd^ G - / fd^ for 

the cubature formula. 

Our research may be considered as a contribution to the topic of ana- 
lytic continuation of solutions to elliptic equations (in particular, harmonic 
functions) , see the discussion and references to the works of V. Avanissian, 
P. Lelong, C. Kiselman, J. Siciak, M. Jarnicki, T. du Cros, on p. 54 — 55 in 
[2], [4], [34], [13], [22], [25]. On the other hand, the concept of polyharmonic 
Hardy spaces which we introduce appears to be a new multivariate concept 
of Hardy space which differs from the existing approaches, cf. [55], [5UJ, |41) . 
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[4"2"] . [8], [H], [35]. This notion will be given a thorough study in a planned 
monograph [53]. Another setting where the polyharmonic Gauss- Jacobi cu- 
bature formulas and Hardy space H\ may be naturally developed, is the strip 
[a, b] x R d , cf. [2T], [24]. 
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